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Abstract. We consider the problem of minimising the n th — eigenvalue 
of the Robin Laplacian in K. . Although for n = 1,2 and a positive 
boundary parameter a it is known that the minimisers do not depend on 
a, we demonstrate numerically that this will not always be the case and 
illustrate how the optimiser will depend on a. We derive a Wolf-Keller 
type result for this problem and show that optimal eigenvalues grow at 
most with n 1//iv , which is in sharp contrast with the Weyl asymptotics 
for a fixed domain. We further show that the gap between consecutive 
eigenvalues does go to zero as n goes to infinity. Numerical results then 
support the conjecture that for each n there exists a positive value of a n 
such that the n th ei genvalue is minimised by n disks for all < a < a n 
and, combined with analytic estimates, that this value is expected to 
grow with n 1//Jv . 
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1. Introduction 

Optimisation of eigenvalues of the Laplace operator is a classic topic in 
spectral theory, going back to the work of Rayleigh at the end of the nine- 
teenth century. The first result of this type is the well-known Rayleigh- 
Faber-Krahn inequality which states that among all Euclidean domains of 
fixed volume the ball minimises the first Dirichlet eigenvalue |24[ [12] \TE[ [19] . 
As a more or less direct consequence of this result, it is possible to ob- 
tain that the second Dirichlet eigenvalue is minimised by two balls of equal 
volume. The case of other boundary conditions has also received some at- 
tention and it has been known since the 1950's that the ball is a maximiser 
for the first nontrivial Neumann eigenvalue [26[ [27] and, more recently, that 
it also minimises the first Robin eigenvalue with a positive boundary param- 
eter [31 [5j [11] , while two equal balls are again the minimiser for the second 
eigenvalue [16]. 

In spite of this, for higher eigenvalues very little is known and even proving 
existence of minimisers in the Dirichlet case poses great difficulties. Bucur 
and Henrot showed the existence of a minimiser for the third eigenvalue 
among quasi-open sets in 2000 [6], and it is only very recently that Bu- 
cur [1] and Mazzoleni and Pratelli [21] proved, independently, the existence 
of bounded minimisers for all eigenvalues in the context of quasi-open sets. 

Moreover, in the planar Dirichlet case and apart from the third and fourth 
eigenvalues, minimisers, are expected to be neither balls nor unions of balls, 
as is already known up to the fifteenth eigenvalue [2]. In fact, it is not to 
be expected that the boundaries of optimisers can explicitly be described in 
terms of known functions either, which means that the type of result that 
one may look for should be of a different nature from the Rayleigh-Faber- 
Krahn type. 

For instance, and always assuming existence of minimisers, there are sev- 
eral qualitative questions which may be raised with respect to this and 
related problems and which include multiplicity issues, symmetry and con- 
nectedness properties, to name just a few. However, even such results might 
not hold in full generality, as some recent numerical results seem to indi- 
cate [2]. 

All of the above issues make this field of research suitable ground for the 
combination of rigorous analytic methods with accurate numerical calcu- 
lations in order to explore the properties of such problems. Indeed, and 
although numerical analysis of eigenvalue problems goes back many years, 
within the last decade there have been several extensive numerical studies 
based on new methods which allow us to obtain insight into the behaviour of 
such problems. To mention just two of the most recent related to eigenvalue 
optimisation, see |23[ [2] for the optimisation of Dirichlet and Dirichlet and 
Neumann eigenvalues, respectively. 

The purpose of this paper is to consider the optimisation of higher eigen- 
values X n of the iV-dimensional Robin eigenvalue problem and analyse some 
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of its properties, combining both the approaches mentioned above. From 
a theoretical perspective, we begin by establishing a Wolf-Keller type re- 
sult, which is needed in the numerical optimisation procedure in order to 
check for non-connected optimal domains. We then consider the asymp- 
totic behaviour of both optimal values of A n , which we shall call A*, and 
the difference between A* +1 and A*. The main result here is the fact that 
A* grows at most with n l ^ N as n goes to infinity, and that the difference 
between optima does go to zero in this limit. Note that this asymptotic 
behaviour for optimal eigenvalues is in sharp contrast with Weyl's law for 
the behaviour of the high frequencies for a fixed domain f2, namely, 



where ujn denotes the volume of the ball of unit radius in M. N and is the 
iV-dimensional volume of fi. Finally, we prove some results regarding the 
behaviour of A n (tQ, a) as a function of the parameters t and a. Although 
intuitively obvious and part of the folklore, their proofs do not seem entirely 
trivial and it is difficult to source them precisely in the literature. Hence we 
have included proofs. 

At the numerical level, our results are obtained using a meshless method 
known as the Method of Fundamental Solutions. Since it is, as far as we 
know, the first time that such a method has been applied to the Robin 
problem, we begin by describing it and stating some basic properties. We 
then present the results of the optimisation procedure. This allows us to 
conclude numerically, as was observed in [171 Sec. 3], that the optimiser will 
depend on the value of a for n larger than two, and provides support for the 
conjecture that for small positive values of a the n th eigenvalue is minimised 
by n identical balls. In fact, and assuming that the domain comprising n 
equal balls stops being a minimiser when its n th eigenvalue becomes larger 
than that of the set formed by n — 3 small balls and a larger ball, we show 
that the value of a at which this happens is increasing with n and grows to 
infinity. 

The paper is divided into three parts. In the first we present the analytic 
results described above, together with the corresponding proofs. This is fol- 
lowed by a description of the numerical method used, and finally we present 
the numerical results obtained. 

2. Theoretical results 

We write the eigenvalue problem as 

— Au = An in fi, 



av 

where v is the outer unit normal to f2 and the boundary parameter a > 
is a constant. We will assume throughout this section that Q C R-^ is a 




as n — > oo 





on dQ 
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bounded, open set with Lipschitz boundary, not necessarily connected, with 
iV-dimensional volume |f2| equal to some fixed constant V > 0. We will also 
use a to denote surface measure. As is standard, we will always interpret the 
problem (JTJ) in the weak sense, so that an eigenvalue A 6 R and associated 
eigenfunction u E solve the equation 



for all v G H 1 ^). It is well known that for each Q C M and a > 0, 
there is a discrete set of eigenvalues {\ n (Q, a)}n>i> an positive, ordered by 
increasing size, and repeated according to their respective multiplicities. For 
each n > 1, we are interested in the quantity 



for each fixed V > and a > 0, where we assume f2 belongs to the class of 
all bounded, open, Lipschitz subsets of R N , as well as the properties of any 
associated minimising domain(s) CI* = Cl*(n,V,a,N). As for the Dirichlet 
problem, when n = 1, the unique minimising domain is a ball p3[TT], w hile 
for n = 2 it is the union of two equal balls |16^ [T7] . Unlike in the Dirichlet 
case, no existence result is known for any n > 3; in R 2 , it was shown in [T7] 
that, for each n > 3, there cannot be a minimiser independent of a > 0. 
As the dependence of A*(V,a) on a > is one of the principal themes of 
this paper, we note here the following basic properties of this function. The 
proof will be deferred until Sec. 12.41 

Proposition 1. Let V > and n > 1 be fixed and for each a > let A* (V, a) 

be given by ([3]). Then as a function of a £ [0, oo), A*(V, a) is continuous 
and strictly monotonically increasing, with A*(V, 0) = and A*(V,a) < 
A*(V, oo), the infimal value for the corresponding Dirichlet problem. 

Remark 2. Throughout this section, we will tend to assume for simplicity, 
especially in the proofs, that (J3j) does in fact possess a minimiser fl*, for each 
n > 1 and a > 0. This assumption can easily be removed by considering 
an arbitrary sequence of domains Qf, with \ n (Q* k ) — > A* . As this type of 
argument is quite standard, we omit the details. Note that without loss of 
generality each domain 0£ may be assumed to have at most n connected 
components, as more could only increase X n (see [TTJ Remark 3.2(h)]). In 
fact, such a sequence may be assumed to be connected. This is a conse- 
quence of results on the stability of solutions to (pQ) with respect to domain 
perturbation: for any domain Q with n connected components, any a > 
and any e > 0, by [TOj, Corollary 3.7], there exists another "dumbbell" -type 
connected domain Q', which has narrow passages joining the disconnected 
components of f2, such that = and A n (fi', a) < A n (f2, a) + e (cf. also 
[TBI Example 2.2]). 

Our point of departure is the way the Robin problem behaves under ho- 
mothetic scaling of the domain. That is, if we denote by tCl the rescaled 




(3) 



A; = A*(V,a) := mf{\ n (Q,a) : \Q\ = V} 
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domain {tx E M : x G O}, then, by a simple change of variables, ([I]) is 
equivalent to 



(4) 



— An = -^-u m 

t 2 

du a „. ^. 

_ + -u = on 9(tfi), 



that is, A n (0,a) = i 2 A n (tf2, a/i), for all (bounded, Lipschitz) Q, C K , 
re > 1, a > and i > 0. (This of course remains valid considering the 
weak form ([2]).) We highlight the change in the boundary parameter. This 

2 

means that, unlike in the Dirichlet and Neumann cases, \tQ,\ n \ n (tQ, a) is 
not invariant with respect to changes in t > 0; rather, the invariant quantity 
is 

(5) \m\ 2 / N x n {tn,a/t). 

This will have a profound effect on the nature of the minimising value A* 
and any corresponding minimising domains. We observe that A*(V, a) as a 
function of a £ (0, oo) may be reformulated as a function taking the form 
A*(iV, a) for t £ (0,oo) and a > fixed, arbitrary. The scaling relation 
gives us immediately that 

(6) t 2 X*(V,-) = \*(r N V,a) 



for all j > 1, all V > 0, all a > and all t > 0, since dU) holds for every 
admissible domain $7 C M. N , so that the same must still be true of their 
infima. Proposition [1] may be reformulated as the following result, which 
will be useful to us in the sequel. The proof will again be left until Sec. 12.41 

Proposition 3. Fix a > andn > 1. As a function ofV 6 (0,oo), AJj(V, a) 
is continuous and strictly monotonically decreasing, with A*(V, a) — > oo as 
V —?■ and X* t (V, a) —?■ as V ->■ oo. 

2.1. A Wolf— Keller type result. An immediate consequence of (|5|) is that 
both the statement and proof of a number of results that are elementary in 
the Dirichlet case now become more involved. Of particular relevance for us 
is the result of Wolf-Keller |28l Theorem 8.1], that any disconnected domain 
minimising A n as in ([3]) must have as its connected components minimisers 
of lower numbered eigenvalues. Here, ([S]) obviously means that we cannot 
hope to be quite as explicit in our description of any potential minimiser. 

Theorem 4. Given V > and a > 0, suppose that there exists a discon- 
nected domain £1* such that = V and A*(V,a) = A n (£l*,a). For every 
1 < k < n — 1, there will exist a unique pair of numbers £1,^2 > 1 (depending 
on k, V, a and N ) with + = 1 which solve the problem 

(7) 

min{max{t?A*(y, t 2 2 X* n _ k (V, £)} :t 1 ,t 2 >l, t^ N + t^ N = 1} 
ti t 2 

= mm{m a x{\* k (q N V,a), \* n „ k (t 2 N V,a)} : t u t 2 > 1, t^ N +q N = l}. 
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Then we may write 

(8) (A£(V,a))S = min \(\* k (V, f ))# + (\* n _ k (V, f))f } . 

Supposing this minimum to be achieved at some j between 1 and n—1, de- 
noting by f^i and Q2 the respective minimisers of X*(V, and A*_ 3 (V, 
/or i/iis pair £2(j) we /iai>e 

(9) n* = U — 2 . 

Because of (0) , we have to define the scaling constants £1 and £2 hi a some- 
what artificial fashion, in terms of a minimax problem (we emphasise that 
Ci)^2 wn l vary with k), and cannot link them directly to the optimal values 
X*j(V,a) as would be the direct equivalent of |28} Theorem 8.1]. Otherwise, 
the proof proceeds essentially as in |28j . 

Proof. We start by proving for each fixed k between 1 and n—1 the exis- 
tence of the pair £1,^2 > 1 as claimed in the theorem. First observe that 
the equivalence of the two minimax problems in ([7]) follows immediately 
from ([6]). Now by Proposition El for V > 0, a > and k > 1 all fixed, as a 
function of t\ £ [1, 00), X k (t^ N V, a) is continuous and strictly monotonically 
increasing from X k (V, a) at t\ = 1 to 00 as t\ — > 00. Moreover, since £2 is 
determined by t\ via the relation ti = (1 — t^[ N )~ l l N , we may also consider 
A* k (t2~ N V,a) as a continuous and strictly monotonically decreasing func- 
tion of t\ G (l,oo], approaching 00 as t\ — > 1 and A*__ fc (V,a) as ii — > 00. 
That is, 

<K-k(0--tI N )~ 1/N v t a) if <i« 1 

> A*_ fc ((1 - t^ N Y l ' N V 1 aj if ti is large enough, 

with the left hand side value strictly increasing and the right hand side 
strictly decreasing in t\. It follows that there exists a unique t\ E (l,oo) 
such that the two are equal. At this value, which we label as t\ =: £1, 
£2 = (1 — ^~ N )~ 1 / N =: £ 2 , the maximum of the two will be minimised. 

Let us now suppose the minimiser Q* of A*(V, a) is a disjoint union 
il* = [/1 U f/2. Since the eigenvalues of £1* are found by collecting and 
ordering the respective eigenvalues of Ui and U2, there exists 1 < k < n — 1 
such that A n (rr,a) = Xk{U\,a). For, if fe = n, then C/2 makes no contribu- 
tion, so rescaling U\ would strictly decrease A n by Lemma [T3l contradicting 
minimality. A similar argument shows that X n (£l*,a) = X n -k(U2,a), since 
otherwise, by expanding U\ and contracting U2, by Lemma [13] we could 
likewise reduce A n (fT,a). It is also clear that Xk(Ui,a) = \* k (\Ui\,a) and 
X n -k(U2, a) = A*_ fc (|C/2|, ct), since otherwise we could replace U\ and/or 
U2 with their respective minimisers and repeat the rescaling argument to 



X* k (t^V,a){ 
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reduce A n (fT,a). Thus we have shown 

X* n {V,a) = X n (n*,a) = X k (U u a) = XU\Ui\,a) = X n . k (U 2 ,a) = X* n _ k (\U 2 \,a) 

We now rescale U\ and U 2 - Let si,s 2 > be such that |siJ7i| = |s2^2| = 
V. Since V = \U\\ + | C/2 1 , we have si,s 2 > 1 and + sj"^ = 1. Now by 



OL . . iV 



(AUV,-))T = (A fc ( fll ?7i,-))T = ( S ^A fc (^i,a))T = ^"(^(V.a))*, 

with an analogous statement for X* n _ k and S2. Adding the two, and using 
that s± N + s% = 1 from the volume constraint, 

f ))* + (Ai_ fc (v, = (A;(y, «))f . 

Sl S 2 

To show that s\ = £1 and S2 = £2, we simply note that, given this A;, the 
unique minimising pair £1,^2 is the only pair of real numbers for which 
£1,6 > 1, Ci N + £2^ = 1 and for which there is equality £i 2 \%(V, §) = 
£ 2 2 A£(V, As si and S2 satisfy exactly the same properties, s\ = £1 and 

Thus we have shown that Q* has the form ([9]), and (jSJ) holds for some 
1 < /c < n — 1 . It remains to prove that AJj is attained by the minimum over 
all such k. To do so, we choose 1 < j < n — 1 arbitrary, label the solution 
to as ji,j2 > 1; and set fi{ to be the domain of volume V such that 

x j (n{,™) = x*(v,-), 
31 J 31 

and analogously for Q? 2 an d A* -(V, a/j±). Now set 

n< = — ni u — ni. 

It is easy to check that = V and that, by choice of ji and j 2 , 

Xj {jKa^j =j?A,- (nj,^) = i 2 2 A, = a, (V 2?a 

meaning that A re (J7j, a) must be equal to all the above quantities. Moreover, 
using the minimising properties of and f^, and that j^[ N + j% N = 1, we 
have 



[X* n (V,a)]- <[X n (n j7 a)] 
proving (JSj) . 



X*[V,- 

31 



+ 



O: 



.72 



□ 
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2.2. Asymptotic behaviour of the optimal values. Another conse- 
quence of © is that any eigenvalue \ n (tft,a) grows more slowly than 
A n ({l) 2 / N as t — > 0. It is thus intuitively reasonable that we might expect 
any optimising domain to have a greater number of connected components 
than its Dirichlet counterpart. Indeed, recalling the variational characteri- 
sation of A n , it is not surprising that increasing the size of the boundary in 
such a fashion carries a fundamentally smaller penalty for the eigenvalues. 
As was noted in [171 Sec. 3], the domain given by the disjoint union of n 
equal balls of volume V/n, call it B n , is likely to play a prominent role in 
the study of A* for sufficiently small positive values of a. Here we go further 
and observe that, simply by estimating Xx(B n ,a), we can already obtain 
quite a strong estimate on the behaviour of A* with respect to n, for any 
a > 0. 

In fact, the following theorem, which again may be seen as an immedi- 
ate consequence of ([5]), shows that we have A* = o(n 1 ^ N+£ ) as n — > oo 
(for any V,a,e > 0), a fundamental divergence from the Weyl asymptotics 
A n (0,a) = 0{n 2 / N ) for any fixed domain O C R N . It is unclear whether 
this is optimal. 

Theorem 5. Given V > and n > 1, let B n denote the domain of volume 
V consisting of n equal balls of radius r = (V/nui^) 1 ^ . Then, for every 
a > 0, 

(10) X* n (V, a) < X n (B n , a) < Na (^) * . 

Proof. Since X n (B n ,a) = Xi(B n ,a), it certainly suffices to estimate the 
latter, that is, to estimate the first eigenvalue of a ball of volume V/n and 
radius r = (V/nu^iy) 1 ^ . Using concavity of Ai with respect to a > 
(Lemma II ip . we estimate this from above by its tangent line at a = (see 
Remark 1 12j) . Since a ball of radius r has volume r N uiN and surface measure 
iVr^-^jv, © at a = gives 

Ai < Ai(fl n ,o)a = 4rrr^ = Nr ~ la - 

| -£»n I 

Substituting the value r = {V /noj^) 1 ^ yields (fTU|) . □ 

2.3. The optimal gap. Adapting an argument of Colbois and El Soufi [7] 
for the Dirichlet case, we may also estimate the dimensionally appropriate 
difference (An+i) N ~ i^n) N ^ 2 f° r each positive V and a, which we do in 
Theorem Such an estimate serves two purposes, giving both a practical 
means to test the plausibility of numerical estimates, and a theoretical bound 
on eigenvalue gaps. In particular, this complements Theorem [5] by showing 
that the optimal gap tends to as n — > oo, albeit not necessarily uniformly 
in a > (see Corollary flOl) . 

The idea is to take the optimising domain $7* for A*, add to it a ball 
B whose first eigenvalue also equals A* and then rescale to obtain a "test 
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domain" for A* +1 . Although the scaling issue ([5]) makes the new behaviour 
possible, it also causes obvious complications, and so we cannot obtain as 
tight a result as in [7j. Instead, we will give two slightly different estimates. 
The first, f)13f) . is tighter but more abstruse, and will be used for computa- 
tional verification, the second one (|14p being more explicit, although only 
in the first case does the bound converge to with n (see Remark [9] and 
Corollary [10]). 

Remark 6. The result that the optimal gap tends to as n — > oo is per- 
haps a priori surprising, and raises the question of whether such a result 
might also be true in the Dirichlet case. Unfortunately, our method tells 
us nothing about the latter, as it rests entirely on the scaling property 
flU). For each fixed n > 1, our bound CH is of the form (X* n+1 ) N/2 - 
{K) N/2 < ^(Ai^i^a))^/ 2 for appropriate constants C = C{N,V) and 
c n = c(N, V, n). The idea is then to show (using that c n — > as n — > oo. 
But fixing n > 1 and letting a — > oo, we recover the bound from [7] of the 
form C(N,V)Xi(Bi), uniform in n > 1. There is no evidence to suggest 
the latter result could be improved. 

The bound (|13p relies on the following auxiliary lemma, needed for a good 
estimate of the constant c n mentioned in Remark [6j 

Lemma 7. Fix V > and a > 0, let A* = A* (V, a) be as in ([3|), and denote 
by B = B(0,r) the ball centred at of radius r > 0. There exists a unique 
value of r > such that 



where ojn is the N -dimensional volume of the unit ball in M and jn_ 

2 

the first zero of the Bessel function Jn__ 1 of the first kind. 



Proof. Consider the left hand side of (|lip as a function of \B\ G (0, oo). An 
increase in \B\ both increases the volume of the domain and decreases the 
Robin parameter. By Lemma \13\ the combined effect must be to decrease 
Ai continuously and strictly monotonically, the latter implying there can be 
at most one value of \B\ giving equality in (jlip . Now note that as \B\ — > 0, 
since V/(V + |£>|) is bounded from below away from zero, \\(B, (V/(V + 



B\))W a ) oo, while if \B\ = V, 

Ai(J3,(l/2)£a) < \i (B, a) < Ai(n*,a) < A n (0*,a) = X* n (V,a), 



where is the minimising domain for ([3]) , and the second inequality follows 
from the Rayleigh-Faber-Krahn inequality for Robin problems [5] Theo- 
rem 1.1]. Hence there must be a value of \B\ in (0, V) for which there is 



(11) 




The corresponding ball B satisfies 



(12) 




2 
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equality in (jlip . To show that the other bound in (112f) also holds, we con- 
sider B r , the ball of radius r = jjv_ 1 1 /-\/A^, where j PiQ is the gth zero of 
the Bessel function J p of the first kind of order p. Then we have 

V 

.V + \B r 

where Xf(B r ) is the first Dirichlet eigenvalue of B r (cf. Lemma ITT]) , and the 
last equality follows from our choice of r. This implies our desired B must 
have radius less than r, giving us (fT2]) , □ 



X\(B r 



a) < Xi(B r ,a) < Af (B r ) = A*, 



Theorem 8. Fix V > and a > and let A* = A£(V, a) be as in ([3]). Lei 

£?* = B*{V,n,a) be the ball satisfying the conclusions of Lemma [7| and let 
B\ denote the ball of unit radius and its N -dimensional volume. Then 



N_ 

2 



(is) (a; +1 ) 

and, weaker but more explicitly, 
(14) 



Ai Si, 



V + \B* 



(\* \t (\*\f^ UJN 



Ai B, 



V \ N 
UJ N 



oo N N a 



a 



Our proof will show, in a manner analogous to the Dirichlet case [7] , that 
given any domain f2 we can find another domain Q of the same volume V 
such that 



X n (£l,a) 



N 

-<^[A l( 5 1)to )]f 



for some appropriate t G (0, 1) inversely proportional to A n (il,a). We omit 
the details. 

Proof. Given V, a > 0, B* and A* as in the statement of the theorem, we 
assume for simplicity that A* is minimised by O*. Let Vt := Q*L)B* (disjoint 
union), set 



V 



V + \B*\ 

and consider the problem ([JJ) on fi, with boundary parameter 

{a on dVl* 

ta on dB*. 

Then by choice of B* and definition of t, 

X n+1 (n,a) = X n (Tt*,a) = Xi(B*,ta). 

Let us now rescale to t£l. If we set 

f a/t on d(tn*) 

a on 9(£B*), 



a :- 
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we have 



Since \t£l\ = V, 

(15) a; +1 < X n +i(tn, a) < \ n+1 (th, a), 

the latter inequality holding since t < 1 so that a < a = a{x) at every point 
of d{tQ). Hence 

1 1 * 1 

K+i < ^A n+ i(0,a) = -p\ n {®> ) = 72 M- 6 >*«)• 

Raising everything to the power of N/2, and subtracting (A* ) N ^ 2 = {\\{B* ,ta)) 



X + i) 2 



(^-l)(X 1 (B*,ta))^ 



Recalling the definition of t, we have 1/t — 1 = |-B*|/y. We will now rescale 
B* , replacing it with a ball of unit radius. That is, letting B* have radius 
r > 0, so that \B*\ = r N uj^, and letting B\ denote the ball of centre and 
radius 1, 



(AiCBVa))T 



Writing 
(16) 



r~ 2 \i{B x ,rta)) 
V\B*\ 



rt 



V + \B*\ 



N 



N 



and \B*\/V = r N uj N /V , we obtain 



(^n+l) 



Ai B 



V + LB* 



which is (|13[) . To remove the explicit dependence on £>* and thus obtain 
(|14p . we simplify the expression (|16p using the crude bounds \B*\ < V in 
the numerator and \B*\ > in the denominator, giving rt < V^^^uj^ 1 ^^ . 
Monotonicity of Ai with respect to the Robin parameter, Lemma [TTT now 
gives CH]). □ 

Remark 9. (i) The bound f)14[) is as explicit as possible for the Robin 
problem, Ai(i?i,/3) being given as the square of the first positive solution of 
the transcendental equation 

^ ^jv(vAi) 



Jn_ 

2 



-i(v / A 1 ) : 



where J p denotes the pth Bessel function of the first kind. 

(ii) In the Dirichlet equivalent of Theorem [SJ the bound is optimal exactly 
for those values of n for which the optimising domain for A* +1 is obtained 
by adding an appropriate ball to the minimiser of A n , which is believed to be 
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true only when n = 1. In our case, since everything converges to its Dirichlet 
counterpart as a — > oo, f)13f) and (|14j) are at least asymptotically sharp for 
n = 1. Moreover, for every n > 1 the bounds converge to zero as a — > 
(as we would hope given that A* — > as a — )■ for every ?i > 1). However, 
even for n = 1, the scaling issue makes it essentially impossible to obtain 
a precise bound for any particular a > 0. Taking N = 2 for simplicity, we 
have A?j = 2A^ in the Dirichlet case, but in the Robin case A?; (a) < 2\\(a) 
for all a > 0, since, denoting by B the ball that minimises AJ(a), we have 
\*(a) = Ai(2~ 1 / 2 5,q) = 2Ai(£,a/2) < 2Ai (£,(*) = 2AJ(a), where we have 
used the fact that the minimiser of A2 (a) is the union of two equal balls [16] , 
the scaling relation @, strict monotonicity of Ai (SI, a) in a (Lemma llip and 
the Rayleigh-Faber-Kralm inequality for Robin problems [5]. Our bound 
in (|13p is, in this case, also smaller than 2\\(a) for all a > 0. However, 
constructing an estimate that involves rescaling domains in this fashion will 
always tend to introduce some error (as happens at (fT5|) in our case), as 
we can never write down explicitly the change in the eigenvalues caused by 
introducing the scaling parameter t into the boundary parameter. 

We now prove the aforementioned result, a complement to Theorem [5l 
that the dimension-normalised gap (A* , i)^ 2 — (A*)^/ 2 approaches zero as n 
goes to 00 for every fixed positive value of V and a. The proof will combine 
(fT3|) with (fT2|) . In the process, we also obtain a growth estimate on A*, but 
this turns out to be weaker than the one found directly in Theorem [5j We 
include the proof of the latter anyway, as both an alternative method and 
to illustrate the principle. 

Corollary 10. For V, a > fixed, as n —> 00 we have 

[A; + i(^»] f -[a;(v»f^o 

and, for every e > 0, 

X* n (V,a) = o(n^ +£ ). 

Proof. Estimating \B*\ from above by (|12|) and from below by in the bound 
(fT3]) gives us 



I 3B—11 

Ai Bi, -2-=^a 



Now, noting that Ai(-Bi, • ) is concave in its second argument, we estimate it 
from above by the corresponding value of its tangent line at a = 0, namely 
T((3) = X' 1 (B 1 ,0)/3. Since Ai(Bi,0) = a(dBi)/uj N (see Remark^, 

TV 



(17) (a; +1 ) t - (K)% < v- 1 ^ 1 -^ (*(dBi)aj»_ ltl ) 2 (a;; 



4 



Observe that, for fixed V, a > 0, the right hand side of (jT7|) converges to 
as n — > 00, proving the first assertion of the corollary. To see the growth 
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bound, label the coefficient of (A*) N ' 4 in (fT7|) as C(N, V, a). Summing 
over n, the left hand side telescopes to give 



(18) {\* n+l )" <(\l)^+C(N,V,a)^2(\ty 



AT 
' 4 



k=l 



Suppose now that for some 7 £ 1, {K) N/2 + 0(n 7 ) as n 00, that is, 
lim suPjj^oq {\* n ) N l 2 /rC 1 = 00. Since A* is increasing in n > 1, a standard 
argument from elementary analysis shows that in fact lim„_ !>00 (A* ) N / 2 /n 1 = 
00, that is, for all Co > 0, there exists uq > 1 such that (A* ) JV//2 > Con 7 for 
all n > no- Hence for Co > fixed, for all n > no we have 



no — 1 n n 

2 

2 



(19) 5>J)-T < M) _T + 0) £ < Ci + C E A: 

fc=l fc=l fc=no fc=l 

where Ci = X^feLT (-^fe) - ^ 4 depends only on N,V,a and no, that is, Co- 
Now we observe 



(20) o(J> 



7/2 \ 

k=l 



0(n 1- *) if 7 e (0,2) 

O(lnn) if 7 = 2, 



(use ELi k ~ S ~ /" if 5 < 1), while ELi fc_7/2 < i+Jr 1 x ~ l/2 dx = 

1 + (7/2 — 1) _1 for all n > 1, if 7 > 2. In particular, combining (|18p and 
(|19p . fixing Co > and a corresponding no > 1 arbitrary, for all n > no, we 
have 



( A* +1 ) - < ( AJ ) - + C(N, V, a) (d + Co ]T k~ * ) , 

k=l 

which we rewrite as 

n 

(21) (A*)f <c 2 + c 3 ^rl 

fe=l 

for all n > no, where C2,C3 > and no > 1 depend only on N, V, a and 
the free choice Co > 0. Recalling that (|2ip holds under the assumption 
(A*)^/ 2 7^ 0(n 7 ) as n — > 00 and using (|2"U|) . this gives us an immediate 
contradiction if 7 > 2/3, forcing (A*)^/ 2 = o(n 2 / 3+e ) for all s > 0. □ 

2.4. Dependence of \ n (tQ,a) on i and a. We will now give some ap- 
pendiceal, but important, results on the behaviour of the Robin eigenvalues 
A n (f2,a) with respect to homothetic changes in $7 or a. Although the ma- 
terial is folklore, we have included a proof as it seems difficult to find one 
explicitly. We will also give the proofs of the corresponding statements for 
A*(V, a), namely Propositions [TJ and 

Lemma 11. For a given bounded, Lipschitz domain $7 C R , and n > 1, 
\ n (ft,a) is an absolutely continuous and strictly increasing function of a £ 
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[0,oo), which is differentiable almost everywhere in (0,oo). Where it exists, 
its derivative is given by 

11 1 1 2 

(22) A An (0,a) = ^^, 

da \\u\\^ n 

where u £ i? 1 ($7) is any eigenfunction associated with A n (Q, a). In addition, 
when n = 1, Xi(Cl,a) is concave, with Xi(Cl,a) < Af(Q), the first Dirichlet 
eigenvalue of CI, and if Cl is connected, then Ai(f2,a) is analytic in a > 0. 

Remark 12. The formula (|22 j) is always valid when n = 1 and a = 0, for 
any bounded, Lipschitz $7 C 1^. In this case, since this corresponds to 
the Neumann problem and the first eigenfunction is always constant, (|22p 
simplifies to 



(23) ^-Ai(n,a) 

da 



_ a(dCl) 

a=0 ~ \Cl\ ' 



a purely geometric property of Cl. The equation (|23|) (which can be obtained 
from a trivial modification of our proof of Lemma [IT]) is reasonably well 
known, and proofs may also be found in \1'6\ I2U|. for example. 

Lemma 13. Given Cl C M. N bounded and Lipschitz, n > 1 and a > 0, 
X n (tCl,a) is a continuous and strictly decreasing function of t £ (0,oo). If 
-jL\ n (Cl, P) exists at ft = ta > 0, then so does 

(24) -X n (m,a) = - ¥ r-j-^-X n (Cl,ta)\, 

where v £ i? 1 (J7) is any eigenfunction associated with X n (Cl,ta). 

Proof of Lemma For the first statement, we note that (weak) monotonic- 
ity and, when n = 1, concavity, are immediate from the minimax formula 
for A n p2 Chapter VI] . We can also derive continuity directly from that for- 
mula, or use the general theory from j!5| Sec. VII. 3, 4]. That is, the form 
associated with (Q]) is 



Qa( n ) = / \Vu\ 2 dx + / au 2 da, 
Jn Jan 



which is analytic in a £ R for each u £ H 1 ^). Hence the associated 
family of self-adjoint operators is holomorphic of type (B) in the sense of 
Kato. It follows that each eigenvalue depends locally analytically on a, with 
only a countable number of "splitting points" , that is, crossings of curves of 
eigenvalues, including the possibility of splits in multiplicities. In our case, 
for each A n (a), the number of such points will certainly be locally finite in a. 
In particular, this means A n (a) consists locally of a finite number of smooth 
curves intersecting each other transversally, so it is absolutely continuous in 
the sense of [251 Ch. 7]. (Throughout this lemma we drop the Cl argument, 
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as it is fixed.) If n = 1 and Q is connected, then \\(a) has multiplicity one 
for all a > and hence no splitting points, so that it is analytic. 

Since [15] also implies that the associated eigenprojections converge, given 
any non-splitting point a (at which A n is analytic), ak — > a and any eigen- 
function u associated with A n (a), we can find eigenfunctions Uk of \ n {a.k) 
such that Uk — > u in L 2 (Q,). We now use a standard argument to show that 
in fact Uk — > u in H 1 ^). Denote by \\v\\* the norm on H 1 ^) given by 
(ll^ u llin + a IMII sn) 1 ^ 2 ' equivalent to the standard one, and assume the 
eigenfunctions are normalised so that ||tt||2,n — H^fclh.Q = 1 for all k > 1. 
Then 



/ \Vu\ 2 + |Vu fc | 2 - 2Vu ■ Vu k dx + / a(u 2 + u\ - 2uu k ) da 
Jn Jan 

= A„(q) + \ n {oLk) - 2A n (o) / uu k dx + (a - a k ) I u\ da, 

Jn Jan 

making repeated use of (|2|). Now we have — > a by assumption, while we 
may use the crude but uniform bound J dn u\ da < \ n {ak)/otk for ak — > a > 
bounded away from zero. Meanwhile, by Holder's inequality, 



uuk dx — I u 2 dx 



as k — > oo, since we know Uk — > u in L 2 (Q,), meaning J n uuk dx — > 1 due to 
our normalisation. Hence, since A n (afc) — > \ n (a) also, 

\\u - life ||* = A n (a) + A n (a fc ) - 2A n (a) / uu k dx + (a-a k ) u 2 k da — >■ 0, 

Jan 

proving Ufe — )• u in both the || . ||*-norm and hence also in the usual ff 1 -norm. 

Let us now compute the derivative of A n (a) at any non-splitting point. 
Suppose a < (3 are two different boundary parameters, with associated 
eigenfunctions u, v G H 1 ^), respectively. Then, using the weak form of A n , 
provided u and v are not orthogonal in L 2 (Q), we get immediately that 

fan uv da 

Xn(j3) - A re a ={p- a) J f - . 

dx 

Now divide through by j3 — a and let /3 — > a. Since we have already seen 
that this forces v — > u in H 1 ^) (also implying that they are not orthogonal 
in L 2 (Q) for f3 close to a), this gives (|22p . In particular, since this is strictly 
positive, and valid except on a countable set of a, we conclude A n is strictly 
increasing. Note that even at splitting points, we may still compute the 
left and right derivatives via this method; we see that it is the change in 
multiplicity (leading to extra eigenfunctions giving different values of (|22p) 
that causes these derivatives to disagree. 

That Ai(n,a) < Af (fi) is an immediate consequence of the minimax 
formulae and the inclusion of the form domains Hq(Q) C if 1 (0). Strict 
inequality is immediate, since Xi(Q,a) is strictly increasing in a > 0. □ 
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Proof of Lemma [TM Since X n (tQ,a) = t~ 2 X n (Q,at), differentiability of the 
former at t is equivalent to differentiability of the latter at (3 = at, and 

^A„(tO, a) = j t (i _2 A n (Q, at)) = -2r 3 A n (ft, at) + a t~ 2 -^X n {Vl, at). 

Using ([22]) and simplifying yields (f!H|) . Continuity of A n (ifi,a) at every 
t > follows immediately from the identity 

X n (sQ, a) — A n (tf2, a) = s _2 A n (0, sa) — i _2 A(f2, ta) 

together with continuity of X n (Q, ta) in t, so that s~ 2 X n (Q, sa) — > t~ 2 X(Q, ta) 
as s — > t. 

Finally, observe that (|24p . holding almost everywhere, also confirms the 
strict monotonicity of \ n (tQ,a) with respect to t > 0. Also note that even 
at points of discontinuity, as in Lemma [TT1 we can again compute left and 
right derivatives, which may disagree due to the change in multiplicity. □ 

Proof of Proposition [7] It follows immediately from the definition of A* (V, a) 
as an infimum and the properties of \ n (£l,a) given in Lemma [11] that 
A*(V, a) is strictly increasing and right continuous in a. Indeed, given 
ao £ [0,oo) and a > ao, if * s a minimising domain so that A*(V, ao) = 
A n (^o,a ), and a > 0, then < A*(V,a) - X n (V,a ) < A*(f2Q,a) - 
X n (Q,Q,ao) — > as a — > ao- For strict monotonicity, if < ao < a, we 
let n* be such that X* n (V,a) = X n (Q*,a) > A n (ft*,a ) > K( V ' a o)- 

Left continuity is harder. We use the property that for each fixed domain 
V, by d22D, 



1 1 2 

(25) — A n (S2,a) = — — < 

eta Il n ll2fi a 

for almost every a > 0. Fixing now ao > and an arbitrary sequence 
< afc < ao, k > 1 with a^ — >■ ao monotonically, we may rewrite (|25p as 

-^-A n (0,a) <CA n (0,a), 
aa 

for almost all a G [ai,ao], where C = 1/ai is independent of Integrating 
this inequality and using the Fundamental Theorem of Calculus applied to 
the absolutely continuous function A n (f2,a) [25] Theorem 7.18], 

X n (n,a ) < A n (ft,a fc )e c(Q °- Qfe) 

for all a/j, k > 1, and all O. Letting fi£ be an optimising domain at a& for 
each k > 1. Then 

A*(V, a ) < liminf A n (f^,a ) 

k— >oo 

<limsu P A n (^,a fc )e ao -^ 

= lim A n (f^,a fc ) = lim A*(V,a fe ). 

k— >oo fc— >oo 

Since this holds for an arbitrary increasing sequence a& — >■ ao, and since the 
reverse inequality is obvious from monotonicity, this proves continuity. 
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Finally, that A*(V, 0) = follows from considering any domain with at 
least n connected components, while to show that A*(V, a) < A*(V, oo), let 
11 be a domain such that \Q\ = V and A*(V, oo) = (Q) > A n (f2, a) > 
A£(V,a). □ 

Proof of Proposition 0. As continuity and monotonicity mirror the proof of 
Proposition [T] closely, we do not go into great detail, but note that now left 
continuity is obvious and the proof of right continuity uses the property (I24p 
to give us the bound 

d 2 

— A n (ifi,a) > -^X n (Q,at) 

in place of ([25]) , If tk — > to is a decreasing sequence, then for C = 2/tf this 
implies 

X n (t n,a) < \ n (t k n,a)e c ^- t0 ^ 
for all k > 1 and all fi, from which right continuity follows in the obvious 
way. 

We now prove that A* (V, a) — > oo as V — > 0. By ©, this is equivalent to 
t 2 A* (1, a/t) — > oo as t ->■ oo. Now 

A;(l,a/0 > Aj(l,a/*) = Ai^, a/t), 

where B' is any ball of volume 1. Since Xi{B',j3) is concave in /3 with 
Ai(S',0) = 0, we have \ x {B',a/t) > Ai(B',a)/t, so that 

t 2 A;(l,a/t) > t 2 Ai(B',a/t) > t\i(B',a) — ► oo 

as t — > oo. Finally, to show that A*(V, a) — > as V — > oo, or equivalently, 
t 2 A* (1, a/t) — > as t — > 0, we simply note that 

t 2 A;(l,a/t) < t 2 A;(l,oo) = C(N,n)t 2 — > 

as t — >• 0, where A*(l,oo) is, as usual, the corresponding infimum for the 
Dirichlet problem. □ 

3. General description of the numerical optimisation 

procedure 

The numerical solution of the shape optimisation problem is divided in 
two steps. At a first level we will describe the application of the Method 
of Fundamental Solutions (MFS) to the calculation of Robin eigenvalues for 
a fixed domain. Then, we will use a steepest descent method (eg. [22]) to 
determine optimal domains for each of the Robin eigenvalues. 

3.1. Numerical calculation of Robin eigenvalues using the MFS. 

We will consider the numerical optimisation of Robin eigenvalues in the class 
of sets S which are built with a finite number of star shaped and bounded 
planar domains. For simplicity, for now, assume that 17 G 3 has only one 
connected component. Thus, f2 is isometric to a domain defined in polar 
coordinates by 

^oo = {(r,0) : < r < roo (t)} 
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where 

oo 

roo(t) = a + (a, cos(i<9) + 6» sin(i0)) . 

8=1 

Now we consider the approximation 

M 

(26) roo(t) ps rjvf (*) := ao + ^ («i cos(i^) + h sin(z0)) , 

i=l 

for a given M E N and the approximated domain £Im denned by 

(27) n M = {(r,9) : < r < r M (t)} . 

Now we describe how to apply the MFS for the calculation of Robin eigen- 
values of fiftf. We take the fundamental solution of the Helmholtz equation 

(28) * A (s) = itf W (VA||s||) 

where |.| denotes the Euclidean norm in IR 2 and Hq is the first Hankel 
function. An eigenfunction solving the eigenvalue problem JT]) is approxi- 
mated by a linear combination 

(29) u{x)^u{x) = Y^^ j (x), 

0=1 

where 

(30) cp J {x) = <5> x {x-y j ) 

are N p point sources centered at some points yj placed on an admissible 
source set which does not intersect Cl (eg. [I]). Each of the point sources <f)j 
satisfies the partial differential equation of the eigenvalue problem and thus, 
by construction, the MFS approximation also satisfies the partial differential 
equation of the problem. We take N p collocation points Xj, i = l,...,N p 
almost equally spaced on 8Qm an d for each of these points we determine 
the outward unitary vector rii which is normal to the boundary at X{. The 
source points j/j are calculated by 

(31) yi = Xi + 7 m, i= 1, N p 

where 7 is a constant (see pQ for details). This choice of collocation and 
source points is illustrated in Figure Q] where the collocation points on the 
boundary and the source points are shown in black and grey, respectively. 
The coefficients in the MFS approximation (|29|) are determined by imposing 
the boundary conditions of the problem at the collocation points, 

du 
ou 

This leads to a linear system 

(32) = t, 



— (xi) + au(xi) = 0, i = 1, N p . 
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Figure 1. The collocation and source points. 

where 

Aij(X) = -Q^i x i) +a<pj(xi), i,j = l,...,N p , 

writing = [fa ... /3 Np \ and = [0 ... 0] 1 . As in Q, we calculate the 
approximations for the eigenvalues by determining the values of A such that 
the system (|32p has non trivial solutions. 



3.2. Numerical shape optimisation. Now we will describe the algorithm 
for the numerical solution of shape optimisation problem associated to the 
Robin eigenvalues. Each vector C := (do, «i ; b±, &2> ^m) defines 

a domain using (|26p and (I27p . The optimisation problems are solved by 
considering each Robin eigenvalue as a function of C, X n (C), and determining 
optimal vectors C. Now we note that we must take into account the area 
constraint in problem ([3]). For each vector C we define a corresponding 
normalised vector C given by 

C 



C 



The domain associated to C has unit area. Now, for each component of the 
vector C, we define an approximation for the derivative of a Robin eigenvalue 
with respect to this component, given simply by a finite difference 



Xn(Vi) -\n[C 



S 



for a small value e, where V% = C + e e% % i = 1,...,2M + 1 and e\ = 
[1 ... 0] T , e 2 = [0 1 ... 0] T , e 3 = [0 1 ... 0] T , .... We then build the 
approximation of the gradient 

dn = [^' n ,l Ki,2 ■•■ Xi,2M+l] 
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which defines the searching direction for the steepest descent method. We 
start by defining Co = C and calculate the next points C n+ \ solving the 
minimisation problem 

Min^ A n \C n - xd n ^j 

using the golden ratio search (eg. PQ). Once we have calculated the optimal 
step length 5 we define 

C n+ i =C n - 5d n , n = 0, 1, 

An alternative procedure would be to apply an optimisation numerical method 
similar to that studied in considering the area constraint. 

3.2.1. Multiple eigenvalues in the optimisation process. As in [2], in the op- 
timisation process we must deal with multiple eigenvalues. For each n, we 

start minimising A n (^C^j and once we obtain 

A n (c) - A n _i (c) < 9, 

for a small value 9 < 1, we modify the function to be minimised and try to 
minimise 

A n (c\ - uj n -i log f\ n (c\ - A n _i (c\\ 
for a sequence of constants uj n -i \ 0. Then, once we obtain 

A n (tf) - A n _! (c) < 6 and A n _i (c) - A n _ 2 (c) < 9, 
we change the function to be minimised to 

A„ (cj - w n _i log (A n (£) - A n _i (cjj - w n _ 2 log (^A n _i ^Cj - A n _ 2 (^C^ 

for suitable choice of constants o; n „i,a; ri _2 \ and continue applying this 
process, adding more eigenvalues to the linear combination defining the func- 
tion to be minimised until we find the optimiser and the multiplicity of the 
corresponding eigenvalue. 

3.2.2. The case of disconnected domains. In the case for which we have a 
set consisting of several connected components, we simply consider a vector 
C defining each component, such that the sum of the areas is equal to one 
and then perform optimisation on these vectors as described above. The ap- 
plication of the MFS in this case is straightforward, considering collocation 
points uniformly distributed on the boundary of each of the components and 
for each of these collocation points, calculate a source point by (|3ip . Note 
that the parametrisation of domains that we considered limits the possible 
shapes to finite unions of star-shaped components. We know that each do- 
main fi£ has at most k connected components (see Remark [2]) and thus, in 
the shape optimisation of A&, there is only a finite number of possibilities for 
building an optimal disconnected domain. These were studied exhaustively, 
working with a fixed number of connected components at each step and 
using the Wolf-Keller type result above (Theorem [4]) to test the numerical 
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results obtained. This process becomes more difficult for higher eigenvalues 
and in that case it might be preferable to use a level set method, as in [23] . 
for instance. 

4. Numerical results 

In this section we will present the main results obtained from our nu- 
merical study on the minimisation of the first seven Robin eigenvalues for 
domains with unit area. We will write B n as shorthand for the domain of 
unit area composed of n equal balls. It is well known that the first and sec- 
ond Robin eigenvalues are minimised by B\ and E>2, respectively. Figure [2] 
shows the evolution of the optimal values of Ai and A2, as a function of the 
Robin parameter a. 




20 40 60 80 100 



Figure 2. Robin optimisers for Ai and A2. 

In Figure [3] we plot the optimal value of A3. We can observe that there are 
two types of optimal domains depending on the value of a. More precisely, 
for a < Q3 ~ 14.51236, the third Robin eigenvalue is minimised by -B3, 
while for a > 03, the ball B\ seems to be the minimiser. In particular, for 
a = «3 uniqueness of the minimiser appears to fail, the optimal value of A3 
being attained by both domains. Note also that in the asymptotic case when 
a —> 00 this result agrees with the conjecture that the ball is the Dirichlet 
minimiser of the third eigenvalue (cf. |28l [T4] ) . 




a 



20 40 60 80 100 

Figure 3. Robin optimisers for A3. 
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While for A3, the Dirichlet minimiser is also the Robin minimiser for 
en > 03, the situation is different for higher eigenvalues. In Figure [51- left 
we plot results for the minimisation of A4. We have marked with a dashed 
line the curve associated with the fourth Robin eigenvalue of the Dirichlet 
minimiser which is conjectured to be the union of two balls whose radii are in 
the ratio \/ jo,i/ji,i, where jo,i and are respectively the first zeros of the 
Bessel functions Jo and J\ (eg. [H]). We can observe, as was to be expected, 
that it is not optimal. The solid curve below it also corresponds to domains 
built with two balls with different areas, but whose optimal ratio of the 
areas depends on the Robin parameter a. It is this family of domains which 
appears to be minimal for larger a. Again we have a value ~ 16.75743 
for which a < implies that the minimiser is -B4. In Figure H] we plot the 
area of the largest ball in the optimal set f2|, as a function of a S [04, 100]. 
We marked with a dashed line the asymptotic case of the Dirichlet optimiser 

for which this quantity is equal to . 2 0.7174. In Figure fright we 

•?0, 



Figure 4. Area of the largest ball of fi|. 

show results for the minimisation of A5. The curve corresponding to the 
Dirichlet minimiser found in [2] is again marked with a dashed line and the 
dotted line below it represents a family of domains of a very similar shape, 
deforming slowly, which appear to be optimisers at their respective values of 
a > a* ~ 40. Then, for 05 < a < a*, where a§ « 18.73537, the minimiser 
is a set composed by a big ball and two small balls of the same area, which 
corresponds to a union of scaled copies of B\ and B2, and where again the 
optimal ratio of these two areas depends on a. For a < 0:5, A5 is minimised 
by B 5 . 

In Figure El- left we plot the results for the minimisation of Xq. The curve 
corresponding to the Dirichlet minimiser is again marked with a dashed 
line, while the Robin optimisers for large a, again close to their Dirichlet 
counterpart, are marked as a sequence of points below it. For some smaller 
values of a it appears that two balls of the same area are the minimiser and 
for a < Q!6 ~ 20.52358, A6 is minimised by Bq. Figure [6}-right shows a zoom 
of the previous figure in the region obtained for a £ [10, 40]. We can observe 




that there is a particular value of a for which the three curves associated 
to unions of balls have an intersection. For this particular a we have three 
distinct optimisers. 

Finally, in Figure [7] we show the results for the minimisation of A7. Again 
the dashed curve is associated to the Dirichlet minimiser and the points 
below it the optimal values obtained for some domains which are very similar 
to the Dirichlet optimiser. 

It is interesting to note that for the region associated with a < 100 
which is plotted in the figure we have a family of domains whose curve is 
below the curve associated to the type of domains which are very similar 
to the Dirichlet optimiser. In the spirit of our Wolf-Keller type theorem 
(Theorem [4]) , this curve corresponds to sets built with an optimiser of \q 
and a small ball, which is an optimiser for Ai. Since the Dirichlet optimiser 



24 P.R.S. ANTUNES, P. FREITAS AND J.B. KENNEDY 




20 40 60 80 



Figure 7. Robin optimisers for A7. 

is connected (cf. [2]), we expect that for some a > 100 the value obtained for 
this type of domain will become larger than the value obtained for domains 
similar to the Dirichlet minimiser, and indeed, testing the algorithm for a = 
300, we find that the optimal domain is connected. For a < a-j ~ 22.167800, 
the minimiser is B7, while for a = aj we have three different types of optimal 
unions of balls. 

4.1. Optimal unions of balls. The above results suggest that the min- 
imiser of A n in two dimensions and for small a consists of the set B n formed 
by n equal balls, and that this situation changes when the line corresponding 
to B n intersects the line corresponding to n — 3 equal balls with the same 
radius and a larger ball, that is, a union of scaled copies of S n _3 and B\. 
The following lemma gives the value at which this intersection takes place 
in terms of the root of an equation involving Bessel functions, showing that, 
for dimension N, this value grows with n 1 ' 1 ^ . 

Lemma 14. Consider the Robin problem for domains in ~ML N of volume V . 
The value a n for which the n th eigenvalue of n equal balls equals the n th 
eigenvalue of the set formed by n — (TV + 1) equal balls with the same radius 
as before and a larger ball is given by 

J f(7o) ,nu> N \i/*f 

2 



OPTIMAL EIGENVALUES OF THE ROBIN LAPLACIAN 25 

Table 1. Optimisers in dimension 2 for a = a n , n = 3, 10. 



n 






3 


14.5123b 




4 


16.75743 




5 


18.73537 




6 


20.52358 




7 


22.16800 




8 


23.69859 




9 


25.13615 


■ o : 


10 


26.49583 





where 70 is the smallest positive solution of the equation 
(33) 

JNA^JpiiCNj) - Cn-jJnJ^Jn.JCnj) + C N 7Jn(j)Jn(C N 7) = 0, 
2 2 2 2 T 22 

with C N = {N + l) l l N . 

Proof. The first eigenvalue of B n , say A, has multiplicity n and is given by 
the smallest positive solution of the equation 

V\Jn (^/\rij — a«/iv_ 1 (v / Ari) = 0, 

where r\ is the radius of each ball, given by {V/{nuj N )yi N . If the smaller 
balls of the domain formed by n — (N + 1) equal balls and a larger ball all 
have radius r\, then the larger ball will have volume (N + l)V/n and radius 
f2 = CVn- We thus have that A equals the n th eigenvalue of this second 
domain, provided that the second eigenvalue of the larger ball, which has 
multiplicity N, also equals A. This is now given by the smallest positive 
solution of the equation 

(1 + ar 2 ) Jn (vAr 2 ) - r 2 v / AJjv. , -, (v / Ar 2 ) = 0. 
2 2 ~ T± 

Writing 7 = y/Xr\, we see that we want to find 7 which is a solution of 

(1 + JXC n1 )Jn{C n1 ) - C n1 Jn +1 (C n1 ) = 

\/~XJn_ (7) — «Jjv_ 1 (7) = 0. 
2 2 

Solving with respect to a/^/\ in the second of these equations, and replacing 
the expression obtained in the first equation yields the desired result. □ 

In the case of dimension two, equation (j33l) reduces to 

JoW-Mv^t) - ^7^0(7)^2(^7) + V3V1 (7)^1(^7) = 
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whose first positive zero is 70 ~ 1.97021, yielding a n « 8.37872- v /ra. Table [Q 
shows the corresponding values of A n for n between 3 and 10, together with 
the different possible unions of balls which also give the same nth eigenvalue 
for a = a n . 

The above numerical results together with Lemma [14] suggest that, at 
least in dimension two, given a fixed value of a there will always exist n* 
sufficiently large such that the minimiser for A„ is B n for all n greater than 
n . 




Figure 8. Numerical values for bound (|13p with n=l,2,...,6. 



4.2. Verification of numerical results using (fl~3j) . As a test of the plau- 
sibility of our numerical results, we finish by computing the error in bound 
(|13p for the case N = 2 and V = 1 . This is shown in Figure [8] where we 
have plotted the quantity 

A* +1 - A* - 7r Ai FBi, ^ j-L 

as a function of a, for n=l,2,...,6, which according to (|13p must always be 
negative. 




5. Discussion 

By combining computational and analytical techniques we were able to 
address the problem of optimizing Robin eigenvalues of the Laplacian, ob- 
taining results for the full frequency range and in general dimensions. The 
application of the MFS to Robin problems provided a fast reliable method 
with which to apply a gradient-type optimization algorithm, yielding min- 
imizers for positive boundary parameters a and up to A7. From this we 
conclude that, except for the first two eigenvalues, optimizers will depend 
on the boundary parameter and approach the Dirichlet optimizer as a goes 
to infinity. 
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In order to address issues related to the connectedness of minimizers we 
derived a Wolf-Keller type result for Robin problems which is also useful to 
identify points where there are multiple optimizers and transitions between 
different branches occur. In particular, as a decreases (while keeping the 
volume fixed), optimizers tend to become disconnected, in contrast with the 
limitting Dirichlet case. From the numerical simulations and the analysis of 
the transition point between the re th eigenvalue of n equal balls and that of 
n — {N + 1) equal balls and a larger ball, we conjecture that for each n there 
exists a transition point, say a n , below which the optimizer for X n consists 
of n equal balls and that a n grows with n 1 ^ . 

Finally, we were able to show that optimizers do not follow Weyl's law 
for the high frequencies, growing at most with n 1 ^ , as opposed to the 
asymptotics for a fixed domain whose leading term is of order r?l N . As far 
as we are aware, it is the first time that such behaviour has been identified. 
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